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Metastatic prostate cancer is currently 
incurable. Metastasis is thought 
to result from changes in the expression 
of specific metastasis-driving genes in 
nonmetastatic prostate cancer tissue, 
leading to a cascade of activated downstream 
genes that set the metastatic process in 
motion. Such genes could potentially serve 
as effective therapeutic targets for improved 
management of the disease. They could 
be identified by comparative analysis of 
gene expression profiles of patient-derived 
metastatic and nonmetastatic prostate 
cancer tissues to pinpoint genes showing 
altered expression, followed by determining 
whether silencing of such genes can lead 
to inhibition of metastatic properties. 
Various hurdles encountered in this 
approach are discussed, including (i) the 
need for clinically relevant, nonmetastatic 
and metastatic prostate cancer tissues 
such as xenografts of patients' prostate 
cancers developed via subrenal capsule 
grafting technology and (ii) limitations in 
the currently available methodology for 
identification of master regulatory genes. 

BACKGROUND 

Metastatic prostate cancer, as distinct from 
localized prostate cancer, is highly resistant 
to conventional therapy. The development 
of new, effective therapeutic strategies 
especially targeting metastatic prostate 
cancer is urgently needed for improved 
disease management.^'^ Metastasis is thought 
to result from changes in the expression of 
specific genes (e.g. upregulation) that lead 
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to activation of multiple downstream genes 
mediating the metastatic process (Figure 1).^'^ 
Such metastasis- driving genes could serve as 
effective therapeutic targets for management 
of metastatic prostate cancer. 

Metastasis-driving genes could be 
identified by comparing gene expression 
profiles of metastatic and nonmetastatic 
human prostate cancer cells in patients' 
specimens to pinpoint genes showing altered 
expression, followed by determining whether 
silencing of such genes can lead to inhibition 
of metastatic properties. A major hurdle 
using this approach, however, is that primary 
prostate cancer samples, the usual source of 
nonmetastatic prostate cancer cells, do not 
consist of pure nonmetastatic cells, but also 
contain metastatic cells. To overcome this 
hurdle, we have developed transplantable. 




^ METASTASIS 



Figure 1: Metastasis is thought to result from 
changes in the expression of specific master 
regulatory genes (e.g. upregulation) that lead to 
activation of cascades of downstream genes that 
set the metastatic process in motion. 



metastatic and nonmetastatic patient- derived 
prostate cancer tissue xenograft lines in 
NOD-SCID (nonobese diabetic/severe 
combined immunodeficient) mice. The 
subrenal capsule grafting technique used to 
develop such cancer tissue lines tends to retain 
important properties of the original cancers, 
including histopathology, chromosomal 
aberrations and gene expression profiles, 
resulting in xenograft lines that closely 
resemble the original cancers.^"^ Recently, we 
have developed paired metastatic (LTL-313H) 
and nonmetastatic (LTL-313B) prostate cancer 
sublines from one patient's specimen, that are 
suitable for investigating the molecular basis 
of metastasis. 

In the present study, comparative analysis 
of gene expression profiles of the LTL-313H 
and LTL-313B sublines led to identification 
of a number of novel, differentially expressed 
genes, including TIMELESS and DLXi, that 
have not previously been associated with 
prostate cancer metastasis and could have a 
metastasis- driving function. 

STUDY DESIGN 

Gene expression profiles of paired metastatic 
LTL-313H and nonmetastatic LTL-313B 
prostate cancer tissue xenografts were obtained 
using microarray technology. Raw gene 
expression data were filtered for improved 
quality prior to analysis of differential gene 
expression. Specifically, probes without 
corresponding gene annotations and 
probes without detectable expression levels 
(less than 3 in log2 scale) were removed. The 
MIAME- compliant gene expression data are 
accessible through GEO Series accession 
number GSE41193 (http://www.ncbi.nlm. 
nih.gov/geo/query/acc.cgi?acc = GSE41193). 

Gene expression profiles of Memorial 
Sloan-Kettering Cancer Center (MSKCC) 
prostate tissues (normal tissues, primary 
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and secondary prostatic adenocarcinomas) 
and associated clinical information^^ were 
downloaded from the CBio Cancer Genomics 
Portal website. Correlations were sought 
between poor prognostic factors of the 
patients (MSKCC cohort) and the relative 
expression levels in their prostate cancer tissues 
of selected genes (identified in the xenografts). 

To validate differential gene 
expressions found, quantitative reverse 
transcription-polymerase chain reaction 
(qRT-PCR) was applied using normal human 
prostate and prostate cancer specimens frozen 
in optimal cutting temperature compounds 
(OCT); RNA was isolated using the RNeasy 
FFPE kit (Qiagen Inc., Hilden, Germany) 
following the manufacturer's instructions. 
The specimens had been obtained from 
patients, with their informed written consent, 
following a protocol approved by a local 
ethics committee (Kelowna Hospital, British 
Columbia, Canada) and examination by a 
pathologist.^^ 

Effects of gene silencing on cell 
proliferation, cell migration and tissue 
invasion were determined using; e.g. PC3M, 
PC3 and C4-2 human prostate cancer cells, as 
previously described. 

Statistical significance was established 
using the Student s t-test 

HIGHLIGHTED RESULTS 

Elevated expression of TIMELESS andDLXl 

Comparative analysis of microarray gene 
expression data (GSE41193), obtained 
from paired metastatic LTL-313H and 
nonmetastatic LTL-313B prostate cancer 
tissue xenografts, was used to identify genes 
in the metastatic tissue that showed elevated 
expression (> 1.5 -fold). For further studies, 
transcription factor-related genes (as defined 
in gene ontology) were selected as they are 
thought to play a critical role in metastasis. ^^'^^ 
Subsequently, genes that had previously been 
associated with metastatic prostate cancer 
were excluded. Of the remaining genes, 
TIMELESS and DLXi, both showing elevated 
expression in the metastatic LTL-313H 
xenografts (vahdated by qRT-PCR; Table 1), 
were of particular interest, since they also 
showed elevated expression in secondary 
prostate cancer samples of the MSKCC cohort 
study. Thus, as shown in Figure 2a, the 
expression level of the TIMELESS gene was 
significantly higher in secondary prostate 
cancer samples in comparison to normal 
prostate and primary prostate cancer samples. 
As well, elevated TIMELESS expression was 
positively associated with poor patients' 
prognostic factors, including increased lymph 



node involvement, increased Gleason score 
and elevated PSA levels (Figure 2b). In the case 
of DLXl , its expression level in the MSKCC 
cohort was significantly higher in both primary 
and secondary cancers when compared to 
normal prostate samples (Figure 2c). This 
differential expression pattern of the DLXl 
gene was validated by qRT-PCR analysis in 
our laboratory using a smaller pool of patients' 
samples (six normal prostate samples and five 
prostate cancer samples; Figure 2d). 

Effect of TIMELESS and DLXl silencing on 
proliferationy migration and tissue invasion 
of prostate cancer cell lines 

siRNA-induced knockdown of TIMELESS 
significantly inhibited the migration 
of PC3M cells in a wound healing 
assay (Figure 3c) without affecting their 
proliferation rate (Figure 3b) or their cell cycle 
progression (data not shown). In contrast. 

Table 1: Elevated gene expression of 
TIMELESS and DLX1 in metastatic LTL-313H 
prostate cancer tissue xenografts compared to 
nonmetastatic LTL-313B counterparts 



Gene 


Fold change 


Fold change 




(microarray data) 


(qRT-PCR validation) 


TIMELESS 


1.83 


1.72 


DLXl 


2.0 


2.89 



DLXl qRT-PCR: quantitative reverse 
transcription-polymerase chain reaction. Fold changes were 
calculated using normalized microarray reads. The data 
were validated by q-RT-PCR 



siRNA-induced knockdown of DLXl did 
not affect the migration or proliferation of 
PC3M (Figure 4) or C4-2 cells (data not 
shown). 

DISCUSSION 

Development of metastasis is generally thought 
to stem from changes in the expression of 
specific, metastasis-driving (master regulatory) 
genes that drive multiple downstream genes to 
set the metastatic process in motion (Figure 1). 
However, the exact molecular and cellular 
mechanisms involved are poorly understood. 
This is in part due to a lack of clinically 
relevant, experimental cancer metastasis 
models. Solid cancers typically consist of 
complex configurations of multiple, diverse 
and interacting cell types. While commonly 
used xenografts of cultured cancer cell lines 
can be useful as models for studying basic 
aspects of the metastatic process, they lack 
the molecular and biological complexity of 
the malignancies (e.g. tumor heterogeneity) 
and hence are not adequate for elucidating the 
multifaceted genetic architecture of metastatic 
cancers and the mechanisms underlying 
their development.^'"^' The LTL-313H and 
LTL-313B prostate cancer xenograft sublines 
used in present study are based on grafting of 
cancer tissues and hence more closely resemble 
the original malignancy. As such, they appear 
to be more suitable for identification of 
prostate cancer metastasis -driving genes. 
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Figure 2: Elevated expression of TIMELESS and DLXl associated with metastatic prostate cancer (PCa). 
(a) Elevated 7"/ME/.ESS expression in metastatic PCa {n = 19) compared to primary PCa (a? = 131) and 
normal prostate tissues {n = 28; microarray gene expression data from the MSKCC Prostate Oncogenome 
Project), (b) Elevated expression of 7"/ME/. ESS associated with poor patient prognosis, (c) Elevated DLXl 
expression in both primary and metastatic tumors compared to normal prostate samples (MSKCC cohort), 
(d) Elevated DLXl expression in five prostate tumors compared to six normal prostate samples. Note that 
three out of six normal samples have undetectable DLXl levels. Statistical significance was established 
using the Student's t-test; *P< 0.05. MSKCC: Memorial Sloan-Kettering Cancer Center. 



Asian Journal of Andrology 



Metastasis-driving genes in prostate cancer 

Chiang Y et al 

547 



While altered expressions of TIMELESS 
and DLXl have been reported for a variety 
of cancers,^""^^ they have not been associated 
with prostate cancer metastasis. The elevated 
expression of TIMELESS found in our 
metastatic xenografts (Table 1), as well 
as in metastatic prostate cancer patients' 
samples (MSKCC cohort; Figure 2a) and 
a strong correlation between elevated 
TIMELESS expression and poor clinical 
patients' outcome (increased seminal vesicle 
invasion and lymph node involvement; 
Figure 2b), suggest that this gene has a role 
in prostate cancer metastasis and perhaps as a 
metastasis- driving gene. However, silencing of 



TIMELESS in PC3M cells only led to limited 
inhibition of cell migration (Figure 3c) and 
had no effect on the metastatic properties of 
C4-2 and PC3 cells, such as cell migration 
and tissue invasiveness (data not shown). 
Similarly, elevated expression of DLXl was 
found to be clinically relevant (Figure 2c 
and 2d), but silencing of DLXl did not show 
any inhibitory effect of metastatic properties 
of the above three prostate cancer cell lines 
in vitro (Figure 4b and 4c). These data 
suggest that TIMELESS and DLXl are not 
metastasis -driving genes. On the other hand, 
the conventional in vitro single cells-based 
assays used might not be adequate (lacking 



biological complexity^^) for evaluating the 
metastasis-driving functions of the genes, 
since it has been reported that metastasis 
can involve the participation of a variety of 
cells, e.g. certain immune cells.^^ In view 
of the above, further studies are needed to 
confirm whether or not TIMELESS or DLXl 
are prostate cancer metastasis -driving genes. 

With regard to alternative approaches 
for identifying metastasis-driving genes, it 
has long been known that master regulatory 
gene networks often act as amplification 
cascades.^^'^^ In such a case, our initial approach, 
aimed at identifying metastasis-driving 
genes on the basis of relatively high gene 
expression changes, would not be effective, 
since changes in the expression of upstream 
master regulatory genes would be much 
smaller than those of their downstream target 
genes. Recently, a promising new approach 
has become available for identification of 
metastasis-driving genes in amplification 
cascades. It is possible to predict upstream 
master regulatory genes through integrative, 
software-based analysis of differential gene 
expression profiles coupled to knowledge of 
upstream regulatory genes (obtained from 
molecular studies). Various programs are 
available such as ARACNE (Algorithm for the 
Reconstruction of Accurate Cellular Networks; 
a novel algorithm, using microarray expression 
profiles) and an Upstream Regulator tool from 
Ingenuity Pathway Analysis.^^"^^ Such strategy 
was effective evidenced by our recent discovery 
of GATA2 gene as a potential metastasis- 
driving gene in prostate cancer. 

Other attempts to identify potential 
master regulatory genes include functional 
genomics approaches. For example, the 
Sleeping Beauty (SB) transposon system 
can be used, in which transposons can 
be randomly inserted into the DNA and 
subsequent development of metastasis can be 
linked to culprit genes.^^ 




Figure 3: siRNA-induced down-regulation of TIMELESS in PC3M cells reduces cell migration. 

(a) Treatment of PC3M cells with siTIMELESS leads to a marked reduction in TIMELESS protein levels, 

(b) but have no effect on cell proliferation, (c) A monolayer of PC3M cells was scratched to examine the 
rate of cell migration into the wounded area. The bar graph represents the percentage of cell-recovered 
wound areas after 24 h of incubation (*P< 0.01). Representative images of the wound captured at 
different time points are shown (at right). Results are representative of three individual experiments 
(mean ± S.D.). S.D., standard deviation. 




Figure 4: siRNA-induced down-regulation of DLXl in PC3M cells has no effect on cell proliferation or migration, (a) Treatment with siDLXl in PC3M cells 
leads to >90% reduction of DLXl gene expression. Downregulated DLXl does not affect (b) cell proliferation and (c) cell migration of PC3M cells. Results 
are representative of three individual experiments (mean ± s.d.). 
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In conclusion, identification of 
metastasis-driving genes is a most challenging 
process. Amongst others, the question is 
raised whether only a transient change in the 
expression of these genes is required to activate 
and maintain the metastatic process, or whether 
their continuous expression is needed. In case 
of the latter, prostate cancer metastasis could be 
arrested by silencing of metastasis-driving genes, 
leading to improved management of the disease. 
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